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ABSTRACT 

We report the first identification of a compact dense galaxy overdensity at z ^ 8 dubbed A2744z80D. 
A2744z80D consists of eight F-dropout galaxies behind Abell 2744 that is originally pinpointed by 
Hubble Frontier Fields studies. However, no studies have, so far, derived basic physical quantities 
of structure formation or made comparisons with theoretical models. We obtain a homogeneous 
sample of dropout galaxies at z ^ 8 from eight field data of Hubble legacy images that are as 
deep as the A2744z80D data. Using the sample, we find that a galaxy surface overdensity value 
of A2744z80D is very high 5 130, where S is defined by an overdensity in a small circle of 6" 

30 physical kpc) radius. Because there is no such a large S value reported for high-z overdensities 
to date, A2744z80D is a system clearly different from those found in previous high-z overdensity 
studies. In the galaxy+structure formation models of Henriques et al. (2015), there exist a very 
similar overdensity, Modelz80D, that is made of eight model dropout galaxies at 2 ; ^ 8 in a 6"-radius 
circle. Modelz80D is a progenitor of a today’s IO^^Mq cluster, and more than a half of the seven 
Modelz80D galaxies are merged into the brightest cluster galaxy (BCG) of the cluster. If Modelz80D 
is a counterpart of A2744z80D, the models suggest that A2744z80D would be a part of a forming 
cluster core of a today’s IO^^Mq cluster that started star formation at 2 ; > 12. 

Keywords: galaxies: formation - galaxies: evolution - galaxies: high-redshift 


1. INTRODUCTION 

Hubble Space Telescope (HST) observations identify 
faint high-redshift galaxies at the epoch of reionization 
and beyond. Hubble Ultra Deep Field 2012 (UDF12) 
project has revealed the faint end of the UV luminosity 
functions of galaxie s at 7 S. z 5, 12 (|Schenker et al.l[2Q13l : 
iMcLnre et al.ll2Ql^ lEllis et al]l2Q13f ). Recently. J. Lotz 
et al. have started the new deep observation program, 
Hubble Frontier Fields (HFF), that carry out deep ACS 
and WFC3-IR imaging in six massive clusters. By us¬ 
ing gravitational lensing magnifications, the HFF project 
finds faint high-redshift galaxies, behind these massive 
clusters, whose intrinsic fluxes reach > 30 magnitude 
(ICoR Rt all lim . Various properties of high-redshift 
galaxies in the HF F data are invest i gated and reported 


by many groups (lAtek et al.l 120151: IZ 

leng et al.l 

12014 

besch et al.ll20l4 llshigaki et al.l 120151: 

Kawamata et al.l 

2015al iMcLeod et al.l 12015 

). Specifically, the studies of 

Ishigaki et al.l (120151) and 

Kawamata et al.l (1201 5aD are 


closely related to this work. 


Three independent studies of HFF have found an over¬ 
density region of 2 : ^ 8 d ropout galaxies in the Abell 
2744 cluster held (A2744C ; Zheng et al.ll2Ql4lAtek et al.l 
I 2 OI 5 I: ITshigaki et al]l2015F Laporte et al.l ( 20141) investi¬ 
gate the HST and Spitzer band photometry of a bright 
source in this region. We refer to this overdensity region 
as A2744z80D. A2744z80D is composed of 8 galaxies 
at 2 : ^ 8 within a circular radius of ^ 6", which cor- 
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responds to ^ 300 comoving kpc in projection. This is 
an overdensity region more dense and compact than the 
previously identified ones at high redshift. For example, 
Brightest-of-Reionizing Galaxies (BoRG) survey reports 
an overdensity of 5 galaxies at 2 : ^ 8 wi thin a mod- 
erately large area of the 127" x 135" sky ([Trent i et al.l 
i^m . Moreover, over densities at 2 : ^ 2 — 7, so far stud¬ 
ied, ar e all defined by large areas with a few como ving 
Mpc^ (iToshikawa et al.ll2Ql4 IVenemans et al.l[200^ see 
Table 5 of iGhiang et al.l I2013D , because these previous 
studies aim to target progenitors of today’s massive clus¬ 
ters. A2744z80D is probably unprecedentedly compact 
and high-dense structure at the redshift frontier of 2 : ^ 8. 
However, no studies have investigated the basic quantita¬ 
tive properties for A2744z80D including galaxy surface 
overdensity. A galaxy surface overdensity 5 is defined by 


where n (n) is the (average) number of galaxies found in 
an area of the sky. In this paper, we evaluate the obser¬ 
vational properties of A2744z80D quantitatively for the 
first time, and discuss the physical origin of A2744z80D 
via the comparisons with the structure formation sim¬ 
ulations of Lambda cold dark matter (AGDM) supple¬ 
mented with galaxy formation models. 

This paper is organized as follows. In Section [2l we 
present details of the observational data. Sample selec¬ 
tion methods are presented in Section [3l and the obser¬ 
vational properties of A2744z80D are described in Sec¬ 
tion IH We compare these properties with the theoretical 
models to identify the model counterpart of A2744z80D, 
and to discuss the physical origin in Section [5l Finally, 
our results are summarized in Section [6l Throughout this 
paper, we adopt a cosmology with = 0.3, Oa = 0.7, 
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and Hq = 70 km s ^ Mpc 

2. DATA 

To estimate 5 of A2744z80D, we use HST deep legacy 
imaging data. We search for useful HST deep legacy 
imaging data, available to date, that are deep enough 
to detect z ^ ^ galaxies as faint as the members of 
A2744z80D with < 29 mag. Among the HST deep 
legacy data, we find that 8 field data sets in HFF, HUDF, 
and their parallel-field regions are useful for our study, 
because these data reach ^ 29 mag and beyond. These 
data sets for our study is presented in Figure [Hand sum¬ 
marized in Table [TJ Note that the data of BoRG and 
CANDELS are too shallow to evaluate 5 of A2744z80D. 
In this section, we explain details of the HFF, HUDF, 
and their parallel-field data sets. 

2.1. HFF 

HFF project plans to observe six cluster and parallel 
fields with ACS and WFC3-IR. However, these obser¬ 
vations are underway. In this study, we use five com¬ 
plete data sets available to date. One is Abell 2744 
cluster (A2744C) data which contains the overdense re¬ 
gion A2744z80D. The others are four parallel field data: 
Abell 2744 parallel (A2744P), MACSJ0416.1-2403 par¬ 
allel (M0416P), MACSJ0717.5+3745 parallel (M0717P), 
and MACSJ1149.5+2223 parallel (M1149P). 

We do not use the associated cluster field data 
of MACSJ0416.1-2403, MACSJ0717.5+3745, and 

MACSJ1149.5-h2223 in the analysis for two reasons. 
The first reason is that the gravitational lensing effects 
distort and narrow the survey area of the cluster fields 
in the source plane. The second reason is that we had 
not made the mass model s of the three clusters yet. 
Recently, iKawamata et al.l (j2Q15b[ ) have provided the 
mass models of the four clusters and the z ^ 8 galaxy 
catalogs in the fields. 

Each data set consists of drizzled science images of 
seven bands in total, three ACS bands, E435W (H 435 ), 
E606W (Veoe), and E814W ( 2814 ), and four WFC3-IR 
bands, E105W (Tios), E125W (J 125 ), E140W (Ji^uo), 
and E160W (i^ieo)- The 5cr limiting magnitudes are 
28.8 — 29.1 in i^ieo band, as summarized in Tabled! 

From September 2013 to February 2014, Frontier Field 
Spitzer program (PI: T. Soifer) has obtained deep IRAC 
images at wavelengths of 3.6 and 4.5/im. This program 
releases corrected Basic Calibrated Data (cBCD) images. 
We use the Spitzer Science Center reduction software 
MOPEX for the reduction of the images. First, we run 
the overlap module, which performs background match¬ 
ing. Then we execute the mosaic module with a drizzled 
factor of 0.65. We obtain the final mosaics with a pixel 
scale of 0".6. The 5cr limiting magnitudes defined with a 
2".0-diameter aperture are 24.7 and 24.9 in the 3.6 and 
4.5 jam bands, respectively. 

2.2. HUDF 

The HUDF is observed by HUDF (PI: Beckwith), 
UDF05 (PI: Stiavelli), UDF09 (PI: Illingworth), and 
UDF12 (PI: Ellis) programs. Combining these data sets, 
the HUDF achieves the quite long exposure time es¬ 
pecially in Tios band, which allows us to make a reli¬ 
able sample of faint galaxies at z ^ 8. HUDF program 


Table 1 

Summary of the HFF, HUDF and the Parallel Data 


Field 

Depth 

( 5(7 limit mag) 

Area^ 

arcmin^ 

Number of 2 : ~ 8 Galaxies 
(same magnitude limit) 

A2744C 

28.90 

1.59^ 

11 

A2744P 

28.75 

3.60 

4 

M0416P 

28.98 

3.81 

8 

M0717P 

28.93 

3.72 

6 

M1149P 

29.07 

3.96 

7 

HUDF 

29.93 

3.74 

7 

HUDFPl 

29.02 

3.72 

5 

HUDFP2 

29.21 

3.04 

7 


^ The depths are defined with a 0735-diameter aperture. 

^ The area of each image is reduced by the conservative masking 
for the areas near the bright foreground stars/galaxies and the 
regions with high noise levels. 

The survey area in the source plane. The area in the image plane 
is 4.79 arcmin^. 

has conducted ultra deep observations with nine bands: 
F775W and F850W in addition to the seven ACS and 
WFC3-IR bands of the HFF survey. We also use the data 
of two HUDF parallel fields, HUDFPl and HUDFP2. 
These fields are observed with the filters F606W, F775W, 
F850W, F105W, F125W, and F160W. The iJieo-band 5cr 
limiting magnitudes are 29.9 in HUDF and 29.0 — 29.2 
in HUDFPl and HUDFP2 (Tabled]). 

3. SAMPLES 

We select z S dropout galaxies based on the Ly¬ 
man break technique. The selection meth o ds ar e ba- 
sical ly the same as thos e of llshigaki et al.l (|2Q15[) (see 
also iHarikane et al.l I2Q15I for the galaxy selection in the 
HUDF, HUDFPL and I IUDF P2 fields). Although Sec¬ 
tion 3 of llshigaki et al.l (j2Q15f ) details the methods, we 
briefly explain them here. 

First, we homogenize the PSF of the images in each 
band in order to accurately compare the magnitudes of 
objects between the different bands. We convolve our 
data with Gaussian kernels to match the FWHMs of stars 
detected in the WFC3 multi-band images. To check the 
reliability of this method, we obtain flux growth curves 
of the stars that are total fluxes as a function of aper¬ 
ture radius in the WFC3 multi-band images. These 
growth curves agree within < 10% in the WFC3 multi¬ 
band images, and the differences of the PSF shapes 
give negligibly small impacts on our color measurements. 
The shapes of PSFs in the ACS images are not impor¬ 
tant, because the ACS images are not used for color 
measurements, but only for flux upper limits, in this 
study. Then we make detection images combining J 125 , 
^^1405 and Hiqq (J 125 and iLieo) images for HFF and 
HUDF (HUDFPl and H UDFP2). We run SExtractor 
(|Bertin fc Arnontill996[ ) in dual-image mode and make 
a source catalog in each field. To measure colors, we use 
MAGJVPER values obtained with a 0735-diameter circular 
aperture. The total magnitudes are estimated from the 
MAGJVPER values and the aperture correction 0.82 mag. 
These total magnitudes agree well with the MAG JVUTO 
magnitudes (see Section 3.1 of llshigaki et al.ll2015f ). 

The selection criteria for z ^ 8 galaxies are represented 
by following equations: 

^105 — Jl25 > 0-5 (2) 

-^125 — ^160 < 0.4. (3) 
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that is the same as the on e used in llshigaki et al.l (|2Q15[ ) 
(see iSchenker et aI]l2Q13l for the original definition). To 
obtain a homogeneous z S galaxy sample from the 
different fields, we apply these selection criteria to the 
source catalogs of all fields, and remove sources whose op¬ 
tical band magnitudes are brighter than the 2cr detection 
levels. We also apply the magnitude limits of J 125 — 29.3 
and Hiqq 29.1 magnitudes in a 0'.'35-diameter circular 
aperture that correspond to the 3.5cr significance levels 
in the shallowest data, and identify 57 galaxies in to¬ 
tal whose positions are indicated in Figure [H We have 
checked the consistency of our catalog with the other 
published ca talogs. Our catalog of A2744C i s cons istent 
with those in I Zheng et al.l (j 2 Ql^ . lAtek et all (j2Q15f ). and 
ICoe et, alJ ll2?fr^. Tn HUDF . we r epro duce sources in 
the ca talog of iMcLure et ^ (|2Q13[ ) and ISchenker et al.l 
(|2Q13[ ) down to the limiting magnitudes of our study. Al¬ 
though the source spatial distribution in M1149P (Figure 
[T|) appears to be biased to the edge of the image, we con¬ 
firm that these are reliable sources in the securely deep 
imaging area by visual inspection. 

Figure [2] presents the surface number densities of these 
galaxies, and indicates that the surface number densities 
are consistent with those of previous studies. The total 
magnitudes are corrected for the magnification by the 
gr avitational lens i ng effe cts. Here we use the mass model 
of llshigaki et al.l (|2Q15l ) which is constructed wi th the 
param etric gravitational lensing package GLAFIC (jOguril 
|2Q1Q[ ). We use magnification-corrected magnitudes for 
the lensed sources in this paper. By these source de¬ 
tection and selection methods applied to all of our data 
sets, we again find the over density of z ^ 8 galaxies, 
A2744z80D, in the north east region of A2744C (Figure 
Pl, as reported by previous studies (| Zheng et al.l l2Ql4 
lAtek et ^l2Q15l : flshigaki et al.ll2Q15[ ). A2744z80D con¬ 
sists of eight z S galaxies within a 6 "-radius circle, 
which a re dubbed ID Yl, Y 2 , Y4, Y5, Y 6 , Y 8 , YIO, and 
Y11 in llshigaki et al.l 1201 51 The magnification factors 
of these galaxies are ^1.5 (jlshigaki et al.l[^Q15F Inter¬ 
estingly, Figure [1] shows no overdense regions similar to 
A2744z80D in the rest of our fields. We confirm that 
A2744z80D is a rare overdensity. 

4. OBSERVATIONAL PROPERTIES OF A2744Z80D 
4.1. Overdensity 

To estimate 5 of A2744z80D, we count the z ^ ^ 
galaxy numbers within a 6 "-radius {'Zi 30 physical kpc) 
circle that corresponds to the angular size of A2744z80D. 
About 16,000 circles are homogeneously placed in our 8 
fields as well as the area centered at A274 4z80D. In the 
A2744 C field, we use the mass model of llshigaki et al.l 
(|2015f ) to correct for the survey area distortion, and mea¬ 
sure the galaxy numbers in the source plane. Figure |4] 
presents the histogram of the galaxy numbers. 

The galaxy-number distribution is basically explained 
by a Poisson distribution shown with the diamonds in 
Figure [H However, the galaxy number in A2744z80D 
which contains eight galaxies significantly exceeds the 
Poisson distribution. If the galaxies are randomly dis¬ 
tributed, the probability to find the A2744z80D excess 
is only 1.9 x 10“^^ in the Poisson statistics. It suggests 
that the galaxy overdensity of A2744z80D is not a sta¬ 
tistical fluctuation, but a physical structure. 


We derive 5 of A2744z80D with Equation ([T]). Follow¬ 
ing the procedures above, we use an area of a 6"-radius 
circle to define galaxy numbers for 5. By these proce¬ 
dures, we obtain 5 = 132^5^ for A2744z80D. Figure [5] 
compares the 5 value of A2744z80D and those of over¬ 
densities at z ^ 2 — 8 reported in previous studies. Here 
we use th e summary of the ove rdens ities at z ^ 2 — 7 pre¬ 
sented in iChiang et al.l (|2Q13f ) and iFranck fc McGaii^ 
( 201511. and the rec ent report of the overdensity at 2 ; ^ 8 
( Trenti et al.ll2Q12l ). ; Because 5 depends on the area used 
for galaxy number counts. Figure [5] shows (5 as a func¬ 
tion of area. Clearly, A2744z80D is placed at a position 
very different from the others (Figure [5]). A2744z80D 
is the most compact and dense high- 2 ; structure, so far, 
identified. 

There are three possible sources of uncertainties in the 
5 estimate. First, there is a possibility that close pair 
sources, Y5-Y8 and Y4-Y6, are not two different galaxies, 
but single galaxies with clumpy structures. The distance 
between Y5 and Y8 is about 0.5 arcsec, corresponding 
to ^ 2.8 kpc projected on the sky, that is the small¬ 
est among the close pairs. The typical size o f 2 ; ^ 8 
galaxies is about 0.5 kpc (jShibuva et al.l [ 20 T 5 I ) that is 
smaller than the Y5-Y8 galaxy distance (^ 2.8 kpc) by 
a factor of 5-6. Some recent studies show that the size 
of hi gh-redshift galaxies at the bright end i s ^1 — 2 
kpc (|Ono et al.ll2013l : iCurtis-Lake et al.ll2016l ). which is 
still smaller than the Y5-Y8 galaxy distance. Thus, it 
is unlikely that the close pairs including Y5-Y8 pair are 
clumps in single galaxies. Second, 2 ; ^ 8 galaxies can be 
blended or obscured by foreground sources, and the de¬ 
tection completeness of 2 ; ^ 8 galaxies may be reduced. 
We calculate the area that foreground sources occupy, 
and find that it is about 4% of the total survey area. 
By the detection incompleteness, the value of 5 would 
change by 4% that is negligibly small. The third source 
of the uncertainty is the foreground source contamina¬ 
tion in our 2 ; ^ 8 galaxy sample. We discuss effects of 
the foreground source contamination in Section 14.21 

4.2. Photometric Redshifts 

To derive the contamination rate of foreground sources, 
we compute the photometric redshifts of our eight 
A2744z80D galaxies. Because rest-frame optical pho¬ 
tometry is key for determining photometric redshifts, we 
measure the source fluxes in IRAC 3.6 and 4.5/im bands. 
Because the IRAC images have a spatial resolution worse 
than the HST images, we model the IRAC images with 
the public softw are for the photome try of low resolution 
images, T-PHOT (jMerlin et al.ll2Q15h . We combine J125, 
^77i4o, and i^ieo images, and use the combined HST im¬ 
age for the high resolution reference image of T-PHOT. 
The source catalog of the reference image is made with 
SExtractor. T-PHOT convolves the reference image with 
the PSE of IRAC, and the convolved reference image is 
fitted to the real IRAC image for each source. We thus 
obtain the best-fit IRAC source fluxes. Eigure [6] shows 
the IRAC magnitudes of the eight A2744z80D galaxies, 
together with the ACS and WEC3 magnitudes. 

We calculate the photometric redshifts with a new ver¬ 
sion of the code H yperz, New-Hyperz, w hich is origi¬ 
nally develo pped bvlBolzonella et al.l (l2QQQl ). We assume 
the S EDs oflBruzual fc Charl(^ (|2QQ3l l models. Chabrier 
IME (jChabrierl l2QQ3f ). exponential SER (r = 0.01 — 3 
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A2744C A2744P M0416P M0717P 



M1149P UDF12 UDFP1 UDFP2 


Figure 1. ifieo images of the 8 fields that we use. The red circles represent the positions of z ~ 8 galaxies. 
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Figure 2. Surface number densites of z ~ 8 galaxies as a function 
of total Hiqo magnitude. The red circles denote the surface number 
densities of our 2 ; ~ 8 galaxy sample that are corrected for the 
leasing effects of the magnification and the distortion in the image 
plane. We show the surface number densities of model 2 ; ~ 8 
galaxies (Section 15.21) with the blue circles, and those of previous 
observational data with the crosses and triangles. 


error of S. 


4.3. Stellar Mass 

We estimate the stellar mass of A2744z80p wit h the 
mass-to-luminosity relation in iLabbe et al.l (j2Q13f ) that 
is derived by S EP fitting of z ^ 8 dropout galaxies. 
Here we assume iChabrierl (j2QQ3f ) IMF in our study to 
obtain stellar masses that can be compar ed with the 
model s sho wn in the next section. Because ILabbe et al.l 
(j2Q13f ) use iSalpeterl (jl955f ) initial mass fun ction (IMF), 
we div ide the mass-to-luminosity relation of iLabbe et al l 
dSoH) by a factor of 1.8: 

M* = 0.9l^;^/1.8 X lO^M© for H - 27.3 (4) 

We confirm that this relat ion is consistent w ith the re¬ 
sults of the recent study (jSong et al.l l2Q15f ). Figure [3 
presents the histogram of the rest-frame UV (i^ieo band) 
total magnitudes for the 8 galaxies in A2744z80D. We 
calculate the total UV luminosity of these 8 galaxies, 
and calculate the total stellar mass with the total UV 
luminosity and the Equation (j4]). The total stellar mass 
is (4 ± 2) X IO^Mq. Note that this total stellar mass is 
comparable to the one of today ’s Large Magellanic Cloud 
(Figure 3 of iDriver et al.l[^Q15f ). We summarize the prop¬ 
erties of A2744z80D in Table [2l 

5. MODEL COMPARISON AND DISCUSSION 


Gyr), metallicity {ZjZ^ = 0.02 — 1.0), Calzetti dust at¬ 
tenuation {Ay = 0.0 — 1.5; iCalzetti et al.l[2000l ). Figure 
[6] present the best-fit SEDs and the redshift probability 
distributions. All the eight A2744z80D galaxies have 
the best-fit SEDs at z ^ 8. Although some faint galaxies 
have moderately high low-redshift contamination rates, 
the total contamination rate is only 16%. This value is 
consistent wit h the e stimate of the contamination rate in 
llshigaki et al.l (|2015f ). This uncertainty is added to the 


5.1. Model 

We investigate theoretical models whether there exist 
a compact and dense overdensity at z ^ 8 similar to 
A2744z80D in galaxy formation models based on the 
ACDM structure formation framework . We use Hen- 
rique s2Q14a galaxy formatio n models (iHenrianes et ^ 
l2Q15f ) made of Millennium (l^ringel et al.l l2QQ5l ) and 
Millennium-II (jBovlan-Kolchin et al.l l2QQ9f ) Simulation 
for the structure formation and the semi-analytic recipes 
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Figure 3. Right: False color image of the A2744C field. The magenta line and circles represent the boundary of the used data and the 
positions of dropout galaxies at 2 ; ~ 8, respectively, in the image plane. The green line and circles are the same, but in the source plane. 
Left: The magnified false color image in the A2744z80D region. The white circles indicate the positions of the 8 ^ ~ 8 galaxies forming 
A2744z80D. 
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Figure 4. Galaxy-number distribution for 2 ; ~ 8 galaxies 
in A2744z80D and the other positions of A2744C, A2744P, 
M0416P, M0717P, HUDF, HUDFPl, and HUDFP2 fields (solid his¬ 
togram). The excess at 8 in the number of galaxies corresponds to 
A2744z80D. The diamonds indicate a Poisson distribution whose 
average number density is the same as our observations. 

for the baryonic processes. Henriques2014a models in¬ 
clude 24 light-cone outputs. Each light-cone has a 3.14 
degree^ area. The light-cone outputs provide catalogs 
with magnitudes of the ACS and WFC3-IR filter system. 

5.2. Model Counterpart of A2744^80D 

Using the catalogs of the light-cones, we select z ^ 8 
dropout galaxies with the color criteria same as our ob¬ 
servations (Section [3]). Figure [2] presents the surface 
number densities of the model dropout galaxies at 2 ; ^ 8. 
The surface number densities of the models are consistent 
with those of observations. The redshift distribution of 
color-selected galaxies in the simulation has the peak at 


00 



4 


Figure 5. Galaxy surface overdensity 5 as a function of area. 
The area is the one used for the S measurements. The red cir¬ 
cle represents A2744z80D at z ~ 8, w hile the blue c ircle d enotes 
the overdensity at 2 ; ~ 8 reported by ITrenti et al.l (120121 b The 
black filled and open circles indicate the o v erden sities at ^ ~ 2 — 7 
summarized in Tab l e 5 o f IChiang et al.l (1201311 and Table 1 of 
IFranck fc McGailghl (120151 b respectively. 


2 ; = 7.8, which is also similar to that in the observation. 
Figure [8] shows the 6 distribution of the model dropout 
galaxies at 2 ; ^ 8, and compares with the one of the obser¬ 
vations, which is the same as Figure |4]but with a different 
scaling for comparison. Figure [8] presents that the S dis¬ 
tributions of the models and the observations are almost 
the same below 4 in number of galaxies. Because the to¬ 
tal area of the observations is significantly smaller than 
the one of the models, the observations do not precisely 
measure the d distribution equal to or greater than 4 in 
number of galaxies, as indicated in Figure IH Thus, Fig¬ 
ures [2] and El suggest that the models well reproduce the 
abundance and clustering of the observed 2 ; ^ 8 galaxies. 

In Figure El we find one region that has the 2 ; ^ 8 
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Figure 6. The SED fits of the eight galaxies in the overdense region. Black arrows present the 2(7 upper limits. We also show the redshift 
probability distributions in each panel. 
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Figure 7. ifieo-band total magnitude histograms of z ~ 8 galax¬ 
ies in A2744z80D (solid line) and the model counterpart Mod- 
elz80D (dashed line). The black filled and open squares indicate 
the mean magnitudes of the A2744z80D and Modelz80D sources, 
respectively. The associated error bars denote the deviation values 
of the magnitude distribution. 
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galaxy number same as A2744z80D. It is located at 
R.A. = 0.64333 deg and Decl. = 0.52677 deg in the cat¬ 
alog Henriques2014a.cones.MRscPlanckl_BC03_010. We 
refer this region as ModelzSOD. Figured presents the UV 
magnitude distribution of ModelzSOD. The mean of the 
UV magnitude distribution of ModelzSOD is similar to 
the one of A2744zSOD within the Icr level. We calculate 
total stellar mass of ModelzSOD from the UV magni¬ 
tudes with Equation |4] in the same manner as Section 
1431 and obtain 2.4 x IO^Mq. This value is comparable 
with the one of A2744zSOD within a factor of 1.5 dif- 



Number of Galaxies 


Figure 8. Galaxy-number distribution of 2 ; ~ 8 galaxies in the 
model catalogs (black histogram). For comparison, we also show 
the galaxy-number distribution from the observations with dashed 
histograms. The dashed histograms are the same as those in Fig¬ 
ure [d] but the total galaxy numbers of the observations are scaled 
to match those of the models. Note that the survey area of the 
simulations is much larger than that of the observations. 

ference. We thus confirm that ModelzSOD is similar to 
A2744zSOD in the UV magnitude distribution and the 
total stellar mass. 

Table [2] summarizes the properties of ModelzSOD. 
Note that there are two total stellar masses for Mod¬ 
elzSOD. One is the estimate with the Equation 0 
(2.4 X IO^Mq), and the other is the value taken from 
the Henriques2014a mod el catalogs (9.2 x IO^Mq) that 
assume the same IME of iChabrierl (|2QQ3f ). Because the 
Henriques2014a model catalog values provide accurate 
stellar masses of ModelzSOD, we hereafter use the stel¬ 
lar mass values in the Henriques2014a model catalogs for 
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ModelzSOD. 

We investigate the structure of ModelzSOD that is se¬ 
lected with a circle of 6"-radms 30 physical kpc) in 
projection. Figure [9] shows the three-demensional map 
in and around ModelzSOD. Seven out of eight galaxies 
of ModelzSOD are physically close to each other in the 
redshift space. The redshift-range width of these seven 
galaxies is Az r\j 0.03 10^ km s ^). This width cor¬ 

responds to ^ 10 comoving Mpc 1 physical Mpc), 
if these galaxies have no peculiar motions against the 
Hubble flow. Because the redshift range of the color se¬ 
lection is A 2 ; ^ 1 (Section [ 3 ]), these seven galaxies are 
very strongly clustered in the redshift space. Figure [9] 
indicates that the distribution of the seven galaxies is 
largely elongated in the line of sight, and that the seven 
galaxies of ModelzSOD are probably found in a filament 
of the large scale structure of the universe. The high 6 
value of ModelzSOD is not only made by the intrinsically- 
high density but by the projection effect. If A2744z80D 
is similar to the model counterpart of ModelzSOD, most 
of the A2744z80D galaxies would form a physical associ¬ 
ation in the redshift space. Moreover, these A2744z80D 
galaxies may reside at a filament of the large scale struc¬ 
tures. 

5.3. Halo and Stellar Mass Evolution 
Suggested by the Model 

We investigate the mass evolutions of the seven galax¬ 
ies of Mode lzSOD that form the physical association 
(Section [Q]) based on Henriques2014a models that trace 
mass evolutions and merger histories of galaxies. Fig¬ 
ure uni presents the halo and stellar mass evolutions of 
the seven ModelzSOD galaxies. Figure flOl indicates that 
the halos of the seven galaxies merge into a single mas¬ 
sive halo with a mass of ^ at z = 0. Because 

this mass is comparable to a galaxy cluster today, Mod¬ 
elzSOD is regard as a cluster progenitor, or protocluster, 
at 2 ; ^ 8. 

The Henriques2014a model catalog has flags indicat¬ 
ing whether a galaxy is included in the central galaxy 
based on the friend-of-friends algorithm. Four out of the 
seven galaxies merge into the central galaxy, the bright¬ 
est cluster galaxy (BCG), that is the most massive galaxy 
in this cluster system at 2 ; = 0. Thus, more than a half 
of ModelzSOD galaxies are building blocks of the BCG. 
In this sense, ModelzSOD is a progenitor of the cluster 
core. Hereafter, we refer to this BGG as ModelzSOD 
BGG. Note that the 2 ; ^ 8 galaxies in the over dense re¬ 
gion do not provide the majority of mass of the BGG, but 
a part of mass of the BGG. At 2 ; ^ 8, the galaxies in the 
overdense region contribute 20% of the total stellar mass 
of the BGG progenitors. The most massive progenitor 
is located about I Mpc away from the overdense region 
at 2 ; ^ 8 (Figure [ 9 ]). Its stellar mass is 2.1 x lO^M©. 
Because ModelzSOD is associated with the most mas¬ 
sive galaxy, it is possible that A2744z80D would be also 
accompanied by such a massive galaxy that would be lo¬ 
cated outside our survey area. This comparison implies 
that such galaxy overdensities can be tracers of massive 
BGG progenitors. 

In order to quantify the relation between overdensities 
and massive BGGs, we calculate the probability that a 
galaxy at 2 ; ^ 8 becomes a BGG in a cluster with the 
halo mass of > lO^^M© at 2 ; = 0. We select model 


-4 -2 0 2 4 2 4 



ARA (cMpc) 


Figure 9. Three-demensional map of model 2 ; ~ 8 galaxies in 
and around ModelzSOD. The filled (open) circles denote positions 
of 2 ; ~ 8 galaxies in (out of) ModelzSOD that is defined by a circle 
of Q" radius on the sky. The white star represents the most massive 
BCG progenitor located 1 Mpc away from the overdense region. 
The top left panel presents the distribution of 2 ; ~ 8 galaxies in 
transverse vs. redshift directions. The bottom panel shows the dis¬ 
tribution projected on the sky. The top right panel is the redshift- 
distribution histograms for the eight galaxies of Modelz80D. 
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Figure 10. Evolutions of the halo masses (dashed lines) and stel¬ 
lar masses (solid lines) evolutions for the seven galaxies in Mod- 
elzSOD. One color of dashed and solid lines corresponds to one 
galaxy. At 2 ; < 4, four out of the seven galaxies merge into the 
most massive progenitor of the BCG. 

galaxies at 2 ; ^ 8 whose UV magnitudes are brighter 
than —20 mag, which is a typical magnitude of galaxies 
in A2744z80D. We find that the probability is 60% for 
the galaxies in overdensities of (5 > 60, and 30% for the 
galaxies in the whole region. The fraction of BCG pro¬ 
genitors in galaxies in overdensities is larger than that in 
the average z ^ S galaxies. 

Finally, we compare the stellar-mass evolution of Mod- 
elz80D BCG with those of the other BCGs in the mod¬ 
els. We select BCGs in clusters whose halo masses are 
^ 10 at z = 0 that is comparable with the one 

of Modelz80D BCG. Figure [TT] presents the stellar-mass 
evolutions of these BGGs. At 2 ; ^ 8, the Modelz80D 
BCG is more massive than the other BCGs at the ^ 2a 
level. Modelz80D BGG is one of the most stellar mas¬ 
sive, i.e. matured, systems among the BGG progenitors 
at 2 ; ^ 8. We investigate the stellar-mass evolution of 
the BGGs to higher and lower redshifts in Figure [TTJ To¬ 
wards low redshift, the stellar masses of the other BGGs 
catch up with the one of Modelz80D BCG at 2 ; < 2. Be¬ 
cause BCGs assemble their masses from a volume larger 
than the 6"-radius circle at z ^ 8, the mass growths 
of some high- 2 ; compact+dense systems such as Mod- 
elz80D could slow down at 2 ; < 2. This indicates that 
the dense systems at 2 ; ^ 8 are not necessarily the pro¬ 
genitors of the most massive systems at 2 ; < 2. Evolution 
to the cluster-sized systems depends on the overdensities 
defined by a circle radius larger than ^ 6" at 2 ; ^ 8. 
Towards high redshift, Modelz80D BGG is significantly 
more massive than the other BGGs at least up to 2 ; ^ 12. 
Although the models have a mass resolution of lO^M© 
that do not allow us to trace back to 2 ; > 12, the mas¬ 
sive system of Modelz80D BGG should start forming at 
a very early epoch. Because S is defined by a very small 
circle of 6" radius (0.3 comoving Mpc) at 2 ; ^ 8 that is 
much smaller than the cluster-scale fluctuations, large 6 
systems like Modelz80D are progenitors of present-day 
massive galaxies that start forming earlier {z > 12) than 
the other massive galaxies. These results suggest a pos¬ 
sibility that overdensities defined by such a small-radius 
circle allow us to identify the site of the earliest star- 


Age of the Universe (Gyr) 



Figure 11. Stellar mass evolutions of BCGs in the Hen- 
riques2014a models. The red line shows the stellar mass evolution 
of the ModelzSOD BCG. Note that ModelzSOD BCG is not asso¬ 
ciated with ModelzSOD at high redshifts. The black dashed line 
denotes the median of stellar masses of the other BCGs which re¬ 
side in clusters with a halo mass of ~ 10^^ Mq at z = 0. The gray 
region represents the 16th-84th percentiles of the BCG stellar-mass 
distribution. 

formation took place in the universe. If A2744z80D is 
truly similar to the model counterpart of Modelz80D, 
the system of A2744z80D would become a cluster to¬ 
day. In this process, some of the A2744z80D member 
galaxies could merge into a BCG to form a core of the 
cluster. We may witness the early cluster core formation 
with A2744z80D. 

These model comparisons provide an interesting in¬ 
sight into future high- 2 ; galaxy over density studies. 
Wide-held observation programs such as Subaru/Hyper 
Suprime-Cam surveys will probably reveal a number 
of compact dense systems at high- 2 ;, wh ich will com¬ 
plement popular studies of protoclnsters (jChiang et al.l 
l2Q13l : iToshikawa et al.l l2Q14l ) dehned by overdensities in 
a large-radius circle. 

6 . SUMMARY 

We investigate the overdense region of eight dropout 
galaxies at 2 ; ^ 8, A2744z80D, behind Abell 2744 clus¬ 
ter, that is originally pinpointed by the several HFF 
studies. This work identihes, for the hrst time, the 
A2744z80D’s physical properties of very compact and 
dense overdensity quantitatively that are clearly differ¬ 
ent from those reported in previous studies of galaxy 
overdensities at 2 : > 2. The detailed comparisons with 
theoretical models are made to understand the physical 
origin of A2744z80D. The results of this study are sum¬ 
marized below. 

1. We make the homogeneous 2 ; ^ 8 galaxy sample 
with the eight field data of Hubble legacy images 
that are deep enough to detect galaxies as faint 
as the members of A2744z80D with < 29 mag, 
and identify A2744z80D in this sample (Figure [3]). 
Moreover, there exist no other over densities sim¬ 
ilar to A2744z80D (Figure [T]). We confirm that 
A2744z80D is a rare overdensity. 

2. We calculate a galaxy surface overdensity S of 
A2744z80D with the homogenous 2 ; r\j 8 galaxy 
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Table 2 

Properties of A2744z80D and ModelzSOD 


5 M. z=o'=’f 

_ [Me] [Mq] [Mq] [Me] 

A2744z80D 8 1321®® (4 ± 2) x 10® 

ModelzSOD 8 159 2.4 x 10® 9.2 x 10® 9.4 x lO^® 3.0 x 10“ 

^ Galaxy number of 2 : ~ 8 galaxies in a 6''-radius circle. 

^ Stellar mass derived from the mass-to-luminosity relation (Equation [4|). 

^ Stellar mass given by the Henriques2014a model catalogs. 

2 : = 0 halo mass taken from the Henriques2014a model catalog. 

^ z = 0 stellar mass taken from the Henriques2014a model catalog. 

^ IChabrieii (|2003) IMF. 


sample. Defining the galaxy surface numbers of 
S with a small circular area of 6" 30 physical 

kpc) radius, we find that S of A2744z80D is very 
large S = 13215^. A2744z80D is the most compact 
and dense galaxy overdensity at z > 2 reported to 
date, and clearly different from those identified in 
the previous studies (Figure [5j). 

3. We estimate the total stellar mass of the eight 
A2744z80D galaxies with the mass-to-l uminosity 
relatio n (Equation|4|) to be 4 x IO^Mq for iChabiiCT 
(|2QQ3[ ) IMF. This total stellar mass is as small as 
the one of today’s Large Magellanic Cloud. 

4. We search for a counterpart of A2744z80D in the 
catalogs of Henriques2014a galaxy formation mod¬ 
els, after we confirm that the models reproduce the 
abundance and the clustering of the observed z ^ 8 
dropout galaxies (Figures [2] and [8]). In the mod¬ 
els, we have found one over density, very similar to 
A2744z80D, that has eight model dropout galax¬ 
ies at z ^ 8 with a similar magnitude distribu¬ 
tion to the one of observations (Figure [7j). This 
model overdensity is referred to as Modelz80D. 
Eight out of seven galaxies in Modelz80D form a 
filament within Az ^ 0.03 (^ 10^ km s“^) that 
is elongated in the line of sight (Eigure [9]). In 
addition to the intrinsically high density in the 
three-dimensional space, this filamentary distribu¬ 
tion makes S of Modelz80D very high. 

5. Investigating the evolution of Modelz80D in the 
models, we have found that Modelz80D is a pro¬ 
genitor of a today’s cluster with a halo mass of 
10 ^"^Mq (Eigure [To]). Moreover, four out of these 
seven Modelz80D galaxies merge into the bright¬ 
est cluster galaxy (BCG) in the today’s cluster. If 
Modelz80D is a similar system to A2744z80D, the 
models suggest that A2744z80D would be a part 
of forming core of the today’s cluster. A2744z80D 
may be a laboratory for the early cluster core for¬ 
mation. Because Modelz80D is a rare (^ 2a) mas¬ 
sive system at least up to 2 ; ^ 12 (Eigure [11]), such a 
compact dense galaxy overdensity would trace the 
galaxy formation in the early universe. 

6. Compact dense high- 2 ; galaxy over densities like 
A2744z80D are unexplored in the past studies 
(Eigure [5]). Because such compact dense galaxy 
overdensities would be clue for understanding early 
cluster core formation as well as the early galaxy 


formation, one needs systematic studies of com¬ 
pact dense galaxy overdensities that will comple¬ 
ment popular studies of protoclusters defined by 
overdensities in a large-radius circle. Euture wide- 
held observation programs such as Subaru/Hyper 
Suprime-Cam surveys will open up the studies of 
compact dense galaxy over densities. 
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